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We investigate the possibility of inducing superconductivity in a graphite layer by electronic corre- 
lation effects. We use a phenomenological microscopic Hamiltonian which includes nearest neighbor 
hopping and an interaction term which explicitly favors nearest neighbor spin-singlets through the 
well-known resonance valence bond (RVB) character of planar organic molecules. Treating this 
Hamiltonian in mean-field theory, allowing for bond-dependent variation of the RVB order param- 
eter, we show that both s- and d-wave superconducting states are possible. The d-wave solution 
belongs to a two-dimensional representation and breaks time reversal symmetry. At zero doping 
there exists a quantum critical point at the dimensionless coupling J/t = 1.91 and the s- and d- 
wave solutions are degenerate for low temperatures. At finite doping the d-wave solution has a 
significantly higher T c than the s-wave solution. By using density functional theory we show that 
the doping induced from sulfur absorption on a graphite layer is enough to cause an electronically 
driven ei-wave superconductivity at graphite-sulfur interfaces. We also discuss applying our results 
to the case of the intercalated graphites as well as the validity of a mean-field approach. 
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I. INTRODUCTION 

The quasi-2D nature of graphite and the linear dis- 
persion with zero density of states a the Fermi level, 
£k = frv/|k|, of a single graphite sheet, called graphene, 
have been known for a long timei More recently, sev- 
eral experimental results have provided evidence that all 
physical properties of graphite are still not yet fully un- 
derstood. For example, traces of superconducting be- 
havior have recently been reported in both graphite 2,3 
and graphite-sulfur composites'^. Highly oriented py- 
rolithic graphite (HOPG) has been shown to display ei- 
ther a partial superconducting or a ferromagnetic-like re- 
sponse to an applied magnetic field even at room tem- 
perature depending on heat treatment and aging^ 3 - For 
graphite-sulfur composites a partially superconducting 
state has been demonstrated below a critical temperature 
as high as T c = 35 K with a continuous cross-over to a 
ferromagnetic-like behavior with increased magnetic field 
or temperature. 3 :^. X-ray analysis of the graphite-sulfur 
composites has shown no significant change in either 
the graphite or the orthorhombic sulfur structure and, 
with the active superconducting volume estimated to be 
only 0.05%, the superconductivity can be understood in 
form of superconducting islands at the graphite-sulfur 
interfaces'^ Large observed magnetic field anisotropy 
strongly suggests that the superconducting correlations 
are localized to the graphite planes in both graphite^ and 
graphite-sulfur composites^. 

In addition to reports on traces of superconductiv- 
ity, both a bulk magnetic-field-driven metal-insulator 
transition (MIT) and a reentrant insulator-metal tran- 
sition (IMT) have been measured in graphite at low 
fieldsi&i Theoretical analysis has suggested that the MIT 



in graphite is due to the magnetic field and/or electronic 
interactions opening up an excitonic gap in the linear 
spectrum, in close analogy with dynamical chiral symme- 
try breaking in relativistic theories of (2+l)-dimensional 
Dirac fermions.— Thermal conductivity data showing 
a kink at the MIT is in agreement with general predic- 
tions of this theoryiii Alternatively, conventional multi- 
band magnetoresistance theory has been used to explain 
the MIT as due to the semimetallic nature of graphite by 
a multiparameter fit to experimental dataJ^ However, 
recent magnetoresistance datai show a non Fermi-liquid 
behavior with respect to magnetic field in the insulat- 
ing phase, thus questioning the applicability of Fermi 
liquid theory. Yet another alternative explanation to 
the MIT in graphite comes from the excellent fifei 3 - to 
the two-parameter scaling analysis suggested by Das and 
Doniach 1 - for a Bose metal-insulator transition. This 
analysis suggests the existence of a non-superconducting 
state of Cooper pairs in the zero temperature limit. In- 
triguingly, the MIT and a suppression of the supercon- 
ducting signal have been found at the same applied field 
in HOPG samples, 3 indicating that superconducting cor- 
relations might be important for the analysis of the MIT. 
Also, the field-induced reentrant IMT has been proposed 
as being caused by Cooper pair formation in the regime 
of Landau level quantization. 6 

With the recent experimental discovery of single layers 
of graphene^ a wide array of theoretical and experimen- 
tal results have also followed on graphene. Although not 
directly connected to what we report here, we mention 
briefly Refs. Il6l|l7l)l8l| l9l | 20l which include measurements 
of proximity induced superconductivity in graphenei 16 i 19 

Our aim in this article is to investigate the effects of 
electronic correlations in a graphite layer. We are espe- 
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cially interested in the possibility to achieve supercon- 
ductivity since many recent experiments seem to point 
to local superconducting correlations in graphite and 
graphitic compounds. We will assume that the inter- 
planar van der Waals coupling in graphite only acts as 
a small perturbation to the in-plane effects and in this 
work exclusively treat only the 2D graphene sheet. The 
carbon atoms in a graphene layer sits in a honeycomb lat- 
tice and are sp 2 hybridized, leaving, in the undoped case, 
one 7r-electron per carbon atom as the valence electrons. 
A well-known model for electronic correlations in the 
half-filled honeycomb lattice is the Hubbard model with 
nearest-neighbor hopping. In the high-[7 limit, the Hub- 
bard Hamiltonian can be rewritten as a t-J Hamiltonian 
where the effective interaction, with coupling constant 
J = At 2 /U, is between nearest-neighbor spins, such as 
to favor resonating valence bond (RVB) spin-singlet cor- 
relations. Due to the large on-site repulsion the Hilbcrt 
space is reduced to not include doubly occupied sites. 
At finite doping, spin-singlet pairs in this model become 
charged superconducting Cooper pairs. Choy et al£^ 
used slave-boson mean-field theory to study the super- 
conducting transition, associated with a non-zero order 
parameter for the RVB spin-singlet, for the t-J model 
in the honeycomb lattice. They found a qualitatively 
similar phase diagram as for the square lattice (sec e.g. 
Ref. I22T ). However, with typical values of the hopping 
t sw 2.5 eV and the on-site repulsion U s» 6 eV for the p z - 
orbitals in graphite, 23 it is questionable to approach the 
effect of correlations in graphite from the high-[/ limit. 
With neither a perturbative weak-coupling approach nor 
the strong-coupling transformation to the t-J model pos- 
sible, theoretical approaches to correlations in graphite 
become highly limited. Early treatments of p7r-bonded 
planar organic molecules and solids such as graphite by 
e.g. Pauling 2 ^ rested heavily on the idea of RVBs, where 
spin-singlet bonds are favored compared to polar configu- 
rations (double and single occupancy of the p z -orbitals). 
For example, good estimates were achieved in this way for 
the C-C bond distance, cohesive energy, and some excited 
state properties. Baskaran^ 5 . used this, half-century old 
concept, in 2002 to propose a phenomenological Hamilto- 
nian for graphite, including nearest neighbor hopping and 
a two-body interaction term favoring nearest neighbor 
spin-singlet formation. The interaction term is identical 
to the RVB term in the t-J model and the only difference 
between the two models is that the phenomenological 
model makes the assumption that double site occupancy 
is still included in the Hilbert space and thus keeps the 
original kinetic energy term. This puts a stronger em- 
phasis on the kinetic term, which in graphite is highly 
non-negligible. But still, the RVB interaction term will 
effectively promote the occurrence of spin-singlet nearest 
neighbor correlations, which, for strong enough coupling 
J (or finite doping), will cause a condensation to a su- 
perconducting state. 

The proposed model is in fact equivalent to the t-J- 
U model with the on-site repulsion U set to zero. This 



model has been suggested as a good candidate for an 
effective one band model of the cuprates^ where d- 
wave superconductivity is triggered by the J interaction 
whereas the U repulsion introduces charge fluctuations 
in the model. It has been shown that a non-zero, repul- 
sive U -term in fact enhances the d-wave superconducting 
correlations on a square lattice i 26 ' 27 i 28 Therefore, under 
the assumption that spin-singlet nearest neighbor cor- 
relations exist in graphite, thus warranting the J-term 
in the Hamiltonian, we expect, by ignoring the on-site 
repulsion U, to use a model which might slightly under- 
estimate the superconducting correlations. 

In this article we will treat Baskaran's effective Hamil- 
tonian in mean-field theory but will allow the order pa- 
rameter to be directionally dependent, allowing for non- 
s-wave superconducting states. We will show that a d- 
wave superconducting state is significantly more favor- 
able than a s-wave state and that only small amounts 
of doping of the graphite layer is enough to initiate the 
onset of superconductivity. We will use this result to ex- 
plore the possibility of <i-wave superconductivity in pure 
graphite and graphite-sulfur composites and also in the 
intercalated graphites (GIC). 

The organization of the paper is as follows. In Section 
II we derive and numerically solve the BCS self-consistent 
equations for the spin-singlet pairing order parameter at 
a specified doping. We deduce all possible symmetries of 
the order parameter as well as briefly discuss their prop- 
erties. In Section III we report band structure calcula- 
tions of graphite-sulfur systems indicating that sulfur can 
induce a hole doping sufficient to cause superconductivity 
in the graphite layer. In Section IV we study the inter- 
calated graphites and show that the spin-singlet <i-wave 
pairing in the 7r-band does not mix with the presumed 
phonon-pairing in the metallic band. We point to dif- 
ferent scenarios which would indicate why the GICs are 
not high-T c d-wave superconductors despite the fact that 
the intercalation induces a high doping into the 7r-band. 
Finally, in Section V we briefly discuss the applicability 
of mean-field theory and summarize our results. 



II. RVB CORRELATIONS IN GRAPHENE 

Motivated by the well-recognized RVB spin-singlet 
character of a graphite layer we model a graphite layer 
using the following Hamiltonian: 25 

H eS = -tJ2(fl9je+9lf j *) 

<i,j>,<7 

i,a <z,j> 

where <i,j> indicates sum over nearest neighbors, fj a is 
the creation operator on the A site of the honeycomb lat- 
tice, and g\ a on the B site, see Figure HJa). The nearest 
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neighbor spin-singlet creation operator is 

4 = ^(44-44) ( 2 ) 

when i £ A-site and has / <-> g for i £ B-site. 

The energy gain from a nearest neighbor spin-singlet 
bond is given by 

-Jhtjhij = J (Si-Sj - 2 n i n JJ ' ( 3 ) 

which clearly displays the similarity between our phe- 
nomenological model and the effective t-J model for 
the high-?/ limit of the Hubbard model. Baskaran 25 
estimated the effective coupling by calculating the 
spin-singlet two-electron ground state energy to J = 
±VU 2 + m 2 - f which with t « 2.5 eV and U w 6 
eV in graphite gives J ft ~ 1. 

To estimate the superconducting transition tempera- 
ture we use mean field theory. A complete Hartree-Fock- 
Bogoliubov factorization will yield order parameters not 
only for the spin-singlet pairing but also for the nearest 
neighbor hopping and site occupation number. Assum- 
ing the latter two to both be bond and site independent, 
we simplify the formalism by assuming that they are al- 
ready included in the effective values of t and /i. Kotliar 
et al£2- showed that for the t-J model on a square lattice 
this is indeed the case even though the superconducting 
order parameter was non-s-wave. With this assumption 
the mean field Hamiltonian can be written in reciprocal 
space as 

H MF = J2 - te4k ' a /iUk CT + H.c. 

k,a,<r 

- £ ^JA Ja e lk a (/t T ^ ki - /t i5 l kT ) + H.c. 

k,a 

+ A^4|A Ja | 2 , (4) 

a 

where Aj a = (/i^j+a) is the order parameter for RVB 
spin-singlet pairing. Here a = ±ai, a2, or a3 denotes the 
three inequivalent directions for nearest neighbors in the 
honeycomb lattice, positive sign for i £ A and negative 
sign for i £ B, see Figure [T] The last constant term is 
irrelevant for the following discussion and will be ignored 
hereafter. The kinetic energy term can be diagonalized 
by the following transformation: 

( ha \ J_ ( rfko- + c kCT \ 

\9kaJ V2 W^(d kff -c kCT ) J ■ {0> 

Here d k creates an electron in the lower 7r-band, c k 
creates an electron in the upper 7r-band, and tpk = 
arg e lk a ) . The fc-dependence of the 7r-bands is given 




FIG. 1: (Color online) Part (a) shows the 2D honeycomb 
lattice with lattice vectors, ci and C2, and nearest neighbor 
vectors, ai, a2, and a3, indicated as well as the two differ- 
ent lattice sites, / (blue) and g (red), (b) Band structure of 
graphene, ±tek, in the first Brillouin zone. The band struc- 
ture close to the corners of the Brillouin zone has a linear 
dispersion with a point-like Fermi surface at zero doping. 

by e k = | Sa e * k a |> see Figure QJb). With this transfor- 
mation the Hamiltonian ^ can be written as 

Huf = E [(/•* + * e k)c k(T c ko - + [fi - te k )d k(T d k(T 

k cr 

- \PlJ A Ja cos(k • a- ^ k )(d kT dL ki - c kT cL ki ) 

a 

- V2i Aja sin(k ■ a - <^ k )(c kT dL ki - d kT c^ ki ) 

a 

(6) 

As seen in Eq. the RVB spin-singlet order parameter 
induces both inter- and intraband pairing in the two it 
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bands. The intraband pairing, second row in Eq. (J6)), 
has an order parameter that is even in fc-space and cor- 
responds to a singlet pairing state. However, the inter- 
band coupling gives a triplet pairing with the so called 
d- vector being d = d z (k)z and with the order parameter 
odd in fc-space. For a bond independent s-wave order 
parameter, the interband order parameter is identically 
zero, but this is not necessarily the case for an arbitrary 
wave symmetry. The Hamiltonian ([6]) can be diagonal- 
ized with a double Bogoliubov-Valatin transformation, 
using two sets of quasi-particle operators: 



Hmf = [ZiAPtoP™ + Aa&O 

k 

+^2,k(7 k0 7ko + 7 k i7ki) 



(7) 



where E^k are the two distinct eigenvalues of Hmf- Close 
to T c the gap parameters Aj a should go to zero and we 
can expand the self-consistency equations in the gap pa- 
rameters to finally arrive at the following self-consistency 
equations at the phase transition temperature: 



— tanh 



2N 
J 



/3 c (te k + fi)\ f/3 c (te^-n) 
1 — tann ' 



k b 



cos(k • a — (fik) cos(k • b — ipjf) 
- sin(k • a — </? k ) sin(k • b — i^k) 



/ tanh(/3 c (£e k + /x)/2) tanh(/3 c (te k - fj)/2) \ 
^ 2(ie k + M ) 2(<e k - A i) J 

sinh(/3 c ^) 

2/icosh(/3 c (ie k + fi)/2) cosh(/3 c (ie k - fi)/2) 

I 



(8) 



'Jb' 



(9) 



Here N is the number of unit cells, (3 C = (fcsT c ) _1 is the 
inverse critical temperature, and 5 is the doping away 
from half-filled 7r-bands in terms of holes per C atom. 
Note however that both self-consistency equations are 
invariant with respect to the sign of the doping, i.e. elec- 
tron and hole doping give the same results. The two first 
terms in Eq. §§§ come from the intraband coupling and 
have the familiar ta,nh(f3E/2)/(2E) BCS-form. The last 
term is due to the triplet interband pairing. For finite 
doping or the reasonably low temperatures we are inter- 
ested in, its contribution is minimal, and we will there- 
fore often ignore this part. For a s-wave solution, i.e. 
Aj a independent of a, Eq. reduces to the compara- 
ble equation found in Ref. l2l|, where the t-J model was 
applied to the honeycomb lattice. For a general solution 
we can, by using O = (A^ ai , A^ a2 , A^ a J T , write Eq. © 
in matrix form as 



J 



-SI 



D B B 
B D B \ 11. 

B B D 



(10) 



where B = B((3 C ) is the RHS of © when a ^ b and 
D = D((3 C ) is the RHS when a = b. The eigenvalues to 
the above equation is easily found to be 



1 
J 



D + 2B, extended s— wave 

D — B, d— wave (and p— wave). 



(11) 



The first solution, the extended s-wave, has eigenvector 
A,/ a = (1, 1, 1) which leaves the intraband order parame- 
ter in Eq. © proportional to e k and the interband order 
parameter identically equal to zero. This extended s- 
wave symmetry belongs naturally to Ai g , the identity ir- 
reducible representation, of the crystal point group Dq^- 
The second solution is two-fold degenerate and spanned 
by vectors A Ja = {(2, -1, -1), (0,-1, 1)}. It leaves both 
the intra- and interband order parameters non-zero. The 
intraband order parameter belongs to the Ei g irreducible 
representation of D§h whereas the interband order pa- 
rameter belongs to the E\ u representation. Both of these 
representations are two-fold degenerate with Ei g repre- 
senting the two four-fold symmetries in a hexagonal crys- 
tal, (k x — fcy, k x ky), i.e. a d-wave state, whereas E\ u has a 
simple odd symmetry, (k Xl k y ), associated with a p-wave 
state. Note that for finite doping or low temperatures, 
the interband pairing is strongly suppressed and the sys- 
tem will then effectively have only a d-wave order pa- 
rameter. Thus, we will hereafter simply call this solution 
d-wave. 



A. Numerical results 

For given i, J, and doping level S (electron or hole dop- 
ing) we can numerically solve Eqs. ([5]) and ©, by sum- 
ming over the irreducible part of the first Brillouin zone, 
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to get the critical temperature T c for both the s-wave 
and d-wave solutions. Figure Ufa) shows the dimension- 
less quantity ksTc/t as a function of J/t for a few dif- 
ferent doping levels whereas Figure dfb) shows T c as a 
function of doping level for t — 2.5 eV and different val- 
ues of J/t around the expected value of 1 for graphenc. 
As seen in Figure for zero doping, the s-wave and d- 
wave solutions are degenerate for low temperatures with 
J/t > J c /t = 1.91 necessary for a superconducting tran- 
sition. This is indicative of a quantum critical point at 
J/t = 1.91 for (5 = 0. This behavior can be understood 
when considering the point-like Fermi surface in graphene 
at zero doping. Our value of J c /t differs from Ref. [25j 
where the critical RVB coupling was estimated using a 
less detailed Cooper pair analysis to J c /t w 3 at zero 
doping. However, with J/t ~ 1 in graphene the qualita- 
tive prediction of no superconducting state in graphene 
is unchanged. 

More important however, is that the d-wave solution 
will always be favored for < 8 < 0.4 because then B in 
Eq. (|10p is negative. In fact, for finite doping the d-wave 
solution has a much higher T c than the s-wave solution. 
While any linear combination of the two d-wave eigenvec- 
tors yields a valid solution at T c , we found numerically 
that only combinations that leave E^ six-fold symmetric 
will be self-consistent solutions below T c . There is only 
one such combination, ft — (1, , ~ 1 ~ v/ ^' ) T , up to 

an arbitrary factor and permutations among the bonds. 
This real-space order parameter leads to complex- valued 
order parameters in reciprocal space and, thus, this so- 
lution is in fact time-reversal symmetry (TRS) break- 
ing. The d-wave solution can be TRS breaking because 
it belongs to a multi-dimensional irreducible representa- 
tion, in this case E2 g , of the D§h group and, if not ig- 
nored, the p-wave solution breaks TRS because it is non- 
unitary. Kuznetsova et al£Q recently pointed out that, in 
the weak-coupling limit, an order parameter belonging to 
Ei g in a hexagonal crystal will be TRS breaking and fer- 
romagnetic. In light of the experimental findings of a pos- 
sible interaction between superconducting and ferromag- 
netic order parameters in graphite-sulfur composites; 3 -^ 
the TRS breaking property of the d- (and p)-wave is very 
intriguing. Note also that this (i-wave solution only has 
a node at T and thus the quasiparticle spectrum is node- 
less. The d-wave solution has the highest T c for 6 = 0.25. 
This is the doping at which the the Fermi level is located 
at the separatrix in the band structure, te^, and thus the 
system has the largest possible Fermi surface. 

As seen in Figure [2Jb) the superconducting state in 
graphene will only be realized for large enough doping 
levels if we assume J/t ~ 1. In recent experiment o 15 ' 31 i 32 
on a single graphene layer doping was induce in the 
graphene by means of applying a gate voltage. However, 
the largest achievable doping level was only 6 ~ 0.008 
electrons/holes per C-atom which is too small for even 
the d-wave solution to have any sizable mean-field su- 
perconducting transition temperature. The experimen- 
tally measured electron and hole concentrations in pure 
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FIG. 2: (Color online) (a) Dimensionless transition tem- 
perature ksTc/t as function of J/t for doping levels S — 
0, 0.001, 0.01, 0.05, 0.1, and 0.2. Big green arrow indicated di- 
rection of increasing doping. The dimensionless temperature 
scale spans approximately 0.03 K to 3000 K for t = 2.5 eV. 
(b) Transition temperature T c calculated using t = 2.5 eV as 
a function of doping level <5 for J/t = 0.8, 0.9, 1, 1.1, and 1.2. 
Big green arrow indicate direction of increasing J/t. 
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graphite, due to c-axis coupling, is even smaller, ~ 10 4 
electrons/holes per C atom. 

III. GRAPHITE-SULFUR COMPOSITES 

The results in the above section indicate that a TRS 
breaking d-wave superconducting state is possible in 
sufficiently electron or hole doped graphene/graphite. 
One such possible system might be the graphite-sulfur 
composites. With the graphite-sulfur composites re- 
taining both graphite and sulfur x-ray signatures and 
the superconductivity possibly located only to the C-S 
interfaces 3 -^ we suggest that it might be a consequence 
of hole doping into the interface graphite layers, which 
then would induce a d-wave superconductivity locally at 
the interface. 

In order to deduce a possible charge transfer between 
C and S we have preformed density functional theory 
(DFT) atomic relaxation and band structure calculations 
of various graphite-sulfur interface systems. Solid sul- 
fur is a molecular solid consisting of Ss rings in an or- 
thorhombic structural unit cell consisting of 16 rings; 33 
Experimentally the sulfur crystal does not have any clean 
surface cuts and we therefore expect the C-S interface to 
be quite disordered. In order to reduce the computa- 
tional cost we have limited our investigation to the in- 
teraction between one graphene layer and various sulfur 
adsorption monolayers. We believe these system to qual- 
itatively capture the graphite-sulfur interface and also 
simulate a realistic experimental situation where atomic 
sulfur is deposited on a single graphene sheet. 

A. Computational details 

In detail, we have studied the following systems: C2S, 
CsS, and CisS with the S atom in various sites on the 
graphene surface. The densest S layer, C2S, has an 
S-S distance of 2.45 A (theoretical lattice constant for 
graphene), in comparison with the 2.08 A distance in 
the Ss ring. The other two systems have the S atoms 
more and more isolated from each other, simulating a 
system of isolated S atoms on a graphite surface. STM 
studies of sulfur electrodeposits on HOPG have shown 
that sulfur is likely to form islands on the graphite sur- 
face in either a C2S configuration or at low concentra- 
tions in a -\/3 x \/3-R30 o lattice corresponding to a CeS 
stoichiometric formula. 34 In addition, any sulfur-graphite 
interface can be expected to contain a fraction of Ss rings 
as well as smaller, reconstructed, rings. We have there- 
fore also studied the interaction of a graphene layer with 
a full Ss ring, with the ring positioned at different sites 
on the graphene surface, in order to capture some effects 
of sulfur molecular interaction with graphite. 

All systems were studied employing the first princi- 
ple density functional theory pseudopotential method 35 
implemented in the Vienna Ab-initio Simulation Pack- 



age (VASP) 3 -^ with the local density approximation 
parametrized by Perdew and Zunger— and projector aug- 
mented wave pseudopotentials 3 ^. The cut-off energy for 
the plane wave expansion was set to 36.7 Ry and wrap- 
around errors were avoided in the FFT-grids. Periodic 
boundary conditions were applied with 22.05 A of vac- 
uum between periodic graphene layers. In all calculations 
the carbon atoms were kept fixed whereas the S atoms 
where allowed to relax in all three spacial directions for 
the C X S systems but in only the z-direction, i.e. perpen- 
dicular to the graphene plane, for the Ss ring systems. 
The atomic structures were relaxed until the forces were 
less than 0.01 eV/A. After k-point convergence tests we 
chose a llxllxl r-centcred grid for the C2S systems and 
7x7x1 grids for the CgS and CigS systems. For the C-Ss 
ring systems we chose a 6x6 graphene unit cell containing 
72 carbon atoms and giving a smallest ring-ring distance 
of 4.05 A. Simulations with only Ss in the same unit cell 
showed that the ring-ring interaction is negligible at this 
distance. We used a 3x3x1 T-centered grid in this case. 



B. Results 

1. C X S 

For the C^S systems we studied three different high 
symmetry positions for S on the graphene surface: Site 
hollow corresponds to S in the middle of the hexagon, 
site top corresponds to S on top of a C atom, and site 
bond corresponds to S on top of a C-C bond. We also 
used several random positions in the xy plane but in no 
case did these sites have a lower energy than the best 
high symmetry site. 

For the C2S system the maximal energy difference be- 
tween the three different sites is very small, 9 meV, with 
the hollow site slightly favored. STM data for the C2S 
system has indicated that S is absorbed in the top or 
bond position^ All systems have a distance between the 
S atom and the graphene surface of 3.3-3.4 A, i.e. much 
larger than the S-S distance, implying that the S atoms 
develop a much weaker bond to the carbon atoms than 
the S-S bond. The CsS and CisS systems have artificially 
much larger S-S distances enforced and as a consequence 
they develop a stronger C-S bond. For both CsS and 
CisS the top positions are favored, only a few meV below 
the bond position but ~ 1 eV below the hollow position. 
In the top and bond cases the S-C surface distance is 1.9- 
2 A. The hollow position retains a larger distance, 2.6- 
2.7 A. These results are in agreement with early quantum 
chemistry cluster calculations, where the top and bond 
positions were favored in several configurations of single, 
double, and triple S atoms configurations^ If we ignore 
interface effects, we can qualitatively compare the ener- 
gies of various systems with the same number of atoms, 
in order to deduce the energetically favorable S densities. 
For example, a system with C2S + 3 x C2 is 2 eV lower 
in energy than CsS, and C2S + 8 x C2 is 2.2 eV lower 
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in energy than CisS. This indicates that sulfur will pre- 
fer to sit on the graphite surface at high densities, and 
if such densities are not available, sulfur island forma- 
tion is likely, as also seen experimentally. 34 We also tried 
to simulate a C2S2 system but this system was found to 
decompose. 

Since at high densities sulfur does not form strong 
bonds to carbon, it is reasonable to assume that the main 
effect of sulfur on the graphene sheet will only be to in- 
duce doping through electron tunneling and not to mod- 
ify the band structure through chemical bond formation. 
This prediction is confirmed in the band structures seen 
in Figure [3] Figure [3^a) shows the band structure along 
high symmetry directions in the first Brillouin zone for 
the C2S system in the hollow structure (red) and pure 
graphene (black). As seen in (a) and the zoom- in around 
the K-point (inset) the sulfur in the C2S system shifts 
the carbon 7r-bands upwards by 0.7 eV and induces only 
a very small hybridization between the carbon and sulfur 
bands. The upward shift leads to an effective hole doping 
in the lower 7r-band, and most importantly, a finite den- 
sity of states at the Fermi level for the lower 7r-band. In 
order to interpret the physical force behind this upward 
shift in energy we also calculated the total energy and 
band structure for a system with the sulfur atoms max- 
imally separated from the graphene layer. In our unit 
cell that leads to a distance of 11 A between the sul- 
fur and graphene sheet. The energy difference between 
this separated structure and the hollow structure proved 
to be very small, 70 meV. Also, the band structure of 
the separated system obviously does not show any band 
mixing but there is still an upward shift of 0.6 eV for the 
7r-bands. Thus we interpret the C2S hollow system as 
a consequence of a physisorption process with weak van 
der Waals bonds between S and C but where the sulfur 
attracts electrons from the graphene sheet through an 
electron tunneling process. The two other systems inves- 
tigated, top and bond, show similar behavior as the hol- 
low structure, but with a somewhat stronger hybridiza- 
tion between the carbon and sulfur. The density of states 
at the Fermi level corresponds in these two cases to an 
upward shift of the 7r-band of about 0.4 eV. This dif- 
ference between the C2S systems can be explained with 
their actual C-S distance, being 3.6 A in the hollow case 
and 3.4 A in the two other cases. 

Figure[3Jb) shows the same results as (a) but for CisS 
in the top (red) position. The top structure has signifi- 
cantly lower energy, 1.4 eV, than the corresponding sep- 
arated structure where C and S are separated by 11 A. 
This is indicative of a chemisorption process with the 
formation of a real chemical bond between carbon and 
sulfur. This strong hybridization creates a 0.3 eV band 
gap in the 7r-band at the Fermi level, and we therefore do 
not have any induced charge transfer in this case. The 
almost energy degenerate bond site also induces a band 
gap for the 7r-band though a bit smaller in size. The band 
structure for CgS is qualitatively the same as for CigS. 

The above energetics and band structures indicate that 




FIG. 3: (Color online) (a) Band structure for C2 (black) and 
C2S in the hollow configuration (red). The Fermi level is 
indicated with a dotted line. Inset shows a zoom-in around 
the K-point. (b) Band structure for Cis (black) and CisS in 
the top configuration (red). Inset shows a zoom- in around 
the F-point. 
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for sulfur atomic absorption on a graphite surface, high 
sulfur concentrations will be preferred, with island forma- 
tion if necessary, and the main effect on the graphite is a 
hole doping of the 7r-band. For the C2S hollow site sys- 
tem the hole doping is 0.015 holes/C-atom. A graphite 
layer with this amount doping will have a transition tem- 
perature to a cZ-wave superconducting state of T c w 10 
K for J/t = 1, but the value is exponentially sensitive 
to the value of the coupling constant </, for example, 
J/t = 1.2 would give T c 80 K. Island formation at a 
lower macroscopic sulfur concentration would mean that 
superconductivity is induced locally in the graphite cov- 
ered by sulfur islands. 



2. C 72 S S 

For the Ss ring systems we studied five different in- 
plane positions allowing the z-direction coordinate, i.e. 
the height of the ring above the graphite layer, to relax. 
As a consistency check we also let one of the five con- 
figurations relax fully, both C and S atoms. The final 
atomic structure showed just a 0.02 A atomic relaxation 
at most and no energy gain. The energy difference be- 
tween all five configurations was small, less than 22 meV, 
with the lowest energies for the rings that had their sul- 
fur atoms far from any underlying carbon atoms. The 
closest distance between the graphene surface and the S 
atoms were in all cases 3.3 — 3.4 A, indicating a weak 
bond between C and S. This could have been expected 
because the sulfur atoms all get a full outer electronic 
shell by the internal bonding in the Ss ring. Since or- 
thorhombic sulfur is a molecular solid we also expect a 
minimal electron tunneling between the graphite sheet 
and the Ss rings. This is confirmed in the band struc- 
ture calculations that show an unaffected and undoped 
carbon 7r-band for all configurations. The ring systems 
are significantly more energetically favorable compared 
to the atomic sulfur deposition configurations, e.g. the 
C72SS system is 6 eV lower in energy than a system with 
8 x C2S + 28 x C2 configuration. This is not very sur- 
prising considering that the Ss ring is the building block 
in solid sulfur and it also means that any intact sulfur 
rings at a graphite-sulfur interface are likely to retain 
their structure. 



IV. INTERCALATED GRAPHITE 

Another graphitic system which has, in fact, been 
known to superconduct for a long time is the intercalated 
graphites (GIC). They are formed by insertion of metal- 
lic atoms in-between the graphite layers. The first dis- 
covered superconducting GICs where alkali metal GICs 
with the composition CsA with A = K, Rb, Cs and with 
T c < 1 K4£ Using pressure during synthesis it is possi- 
ble to increase the metallic concentration and in general 
T c is found to increase with metal concentration up to 



T c = 5 K for C2Na. 41 Recently, the intercalated com- 
pounds CeCa and CgYb were found to be superconduct- 
ing with T c — 11.5 K and 6.5 K, respectively^ 2 . These 
are the highest transition temperatures measured for the 
GICs as of today. 

Intercalating graphite does not significantly change the 
graphite structure but the intercalated species always do- 
nate electrons into the graphite 7r-band so the Fermi sur- 
face of a GIC can schematically be thought of as 7r-band 
cylinders along HKH and a metallic band, often at the 
zone centeri 43 ' 44 However, zone folding will often cause 
the Fermi surfaces to overlap. For simplicity we will here- 
after call the metal band a s-band since that is the case 
in the alkali GICs. The charge transfer from metal to 
graphite is estimated to be be ~ 0.6 for CsK— and re- 
lated compounds and was recently calculated to be 0.32 
for CeCa^. These values give 5 = 0.08 and 0.05 respec- 
tively, which, if only considering the 7r-band, would yield 
a T c of several hundred K for a d-wave superconducting 
state. Not only are these values a magnitude or more 
too large but it has been also found that there is no sim- 
ple relationship between doping and T c in the GICs4£ 
However, in most GICs the metal band is only partially 
occupied and should also be taken into account. While 
exotic pairing mechanisms have been suggested for espe- 
cially the high-T c GICs 47 - several DFT calculations have 
recently found large enough electron-phonon couplings 
to induce the correct T c in C 6 Ca and C 6 Yb4&2£ The 
phonons are found to mainly couple to the s-band since 
the 7r-band will, due to symmetry, not couple to any 
out-of-plane vibrations, therefore limiting the electron- 
phonon coupling for the 7r-band. It is therefore of interest 
to consider the possibility to couple the ci-wave state in 
the 7r-band to a phonon induced s-wave superconducting 
state in the s-band. 

We have investigated the possibility to couple the 
two different order parameter states with both band hy- 
bridization and interband phonons. Both of these ef- 
fects are estimated to be small but we are interested 
in the fundamental possibility to couple the d- and s- 
waves to each other, and possibly suppress the d-wave 
state this way. We consider a mean-field Hamiltonian 
consisting of only the upper 7r-band, since the doping 
is sufficient to effectively suppress the effect of the lower 
band and inter-7r-band coupling, and a BCS-like phonon- 
driven superconducting s-band with order parameter 
A s = X)k( s -k| s kT)- ln the case °f band hybridiza- 
tion, we introduce the additional term (r^c^ Sk CT + H.c.) 
to the Hamiltonian in order to produce a band mixing. 
By re-diagonalizing the kinetic energy, i.e. calculating the 
new band structure with the band mixing included, the 
band operators will change according to c = ac + /3s, 
with a ~ 1 and (3 oc T 2 , and similarly for s. Expressed 
in these new operators the spin-singlet pairing in the 
original 7r-band will induce a s-tt interband pairing as 
well as an additional pairing mechanism for the s-band. 
The final Hamiltonian can be diagonalized with a double 
Bogoliubov-Valatin transformation similar to the proce- 
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dure followed in Section II. The resulting self-consistency 
equations for Aj a and A s have four different solutions of 
which two are the <i-wave states for Aj a with A s = and 
the final two mixes the s-wave state of the spin-singlet 
pairing with the phonon s-wave state. With no coupling 
between the ci-wave state and the s-wave phonon state 
and considering only small values of T we conclude that 
band hybridization, treated in mean-field theory, will not 
alter the fact that the d-wave solution has a significantly 
higher T c than the (coupled) s-wave solutions at the con- 
sidered doping levels. In the case of interband phonon 
pairing we used the two-band model introduced by Suhl 
et al^ but here with a ^-dependent coupling in the 7r- 
band. Even here, the <i-wave solution does not mix with 
the s-wave phonon pairing. In both cases, the lack of 
coupling between the two different order parameters can 
be traced back to the fact that for spherical or 6-fold 
symmetric bands, the two-fold degenerate cZ-wave solu- 
tion will have a zero overlap with any functions whose 
/c-dependence comes solely from the band structures. In 
fact, in order to achieve a non-zero A s even below T c (d- 
wave) the six-fold symmetry of the quasi-particle energy 
has to be broken in order to create non-zero overlap in- 
tegrals. 

The above results indicate that the two-fold degener- 
ate cZ-wave superconducting state in the 7r-band will not 
couple to a phonon-driven superconducting state in the 
metal s-band. With the doping levels relevant for the 
GICs this means that the unaffected mean-field <i-wave 
solution does not agree with experiments. We can see at 
least two possible explanations to this. Quite fundamen- 
tally, with a partially full s-band present, the electronic 
screening of the on-site repulsion in the 7r-band will be 
significantly enhanced. This should mean that the effec- 
tive coupling J is drastically reduced and thus no or only 
a very small effective attraction exists between electrons 
in the 7r-band. Alternatively, or in addition, there might 
be important many-body effects that our in the mean- 
field solution has ignored. Recent band structure calcu- 
lations of GICs have shown that the Fermi surfaces of the 
7r-band and the s-band is likely to overlapj 46 ' 47 One way 
to model the effects of such overlap would be to introduce 
an interband tunneling term of the form Tksl _Ck CT &s m 
the periodic Anderson impurity model but not treat it 
as the band mixing term above. It is important to note 
that such a model for graphite will include a significant 
kinetic energy term compared to a periodic Anderson 
impurity model where the impurity lattice is supposed 
to not have any intrinsic hopping. Colloquially speaking, 
the coupling to an external band will cause fluctuations of 
electrons between the two bands which effectively might 
suppress the spin-singlet interaction in the 7r-band. Such 
an effect might suppress the superconducting correlations 
in the 7r-band to the degree to favor the phonon-driven 
superconductivity in the GICs instead. It is worth not- 
ing that while the atomic sulfur layer in the C2S systems 
show metallic behavior we do not expect a bulk s-band 
since sulfur is a bulk insulator and thus this effect should 



not be present in the graphite-sulfur composites. 



V. CONCLUDING REMARKS 

Assuming there are pronounced spin-singlet correla- 
tions in graphite we have shown that a d-wave TRS 
breaking superconducting state is possible even for only 
moderately doped graphite sheets. The intrinsic carrier 
concentration in graphite due to the interlayer coupling 
as well as reported electrostatically induced doping in sin- 
gle graphene sheets are however both too small to give 
a measurable T c . So far we have, however, completely 
ignored the effect of defects and natural adsorption of 
gaseous molecules in graphitic systems. It has recently 
been shown that in graphene extended defects such as 
dislocations, disclinations, edges, and micro-cracks will 
lead to self-doping where electronic charge is transferred 
to or from the defect to the bulk. 50 Also, natural chemical 
doping has been detected on graphene J£ We expect both 
mechanisms to be more present in graphite samples^ than 
in graphene where the play only a minor role. It should 
then, in principle, be possible to achieve local areas in a 
graphite sample where the electron or hole doping is large 
enough to either cause a condensation of the RVB pairs 
or at least enhance their correlations enough in order to 
measure a superconductor-like response of the sample. 
It is natural that such mechanisms are very sensitive to 
heat treatment as well as aging, as found experimentally. 

In the graphite-sulfur composites the situation is 
slightly more favorable since atomic sulfur will tend to 
hole dope graphite. As we have shown it is possible to 
achieve doping levels of the order of 0.015 holes per C- 
atom in the 7r-band for atomic sulfur deposition on a 
graphite sheet. This is large enough to induce d-wave 
superconductivity with experimentally measurable tran- 
sition temperatures within the graphite layer. In the 
graphite-sulfur composites we anticipate the interface 
to be disordered with occurrence of both atomic sulfur 
and sulfur rings. The sulfur rings do not influence the 
graphite layer to any large extent but the atomic sulfur 
will cause an electron transfer to the sulfur atoms and 
thus induce superconductivity. The anticipated disor- 
deredness of the graphite-sulfur interfaces and the vari- 
able sizes of the graphite and sulfur domains will make 
also these systems sensitive to heat treatment and aging 
as recently demonstrated.— 

We have so far interpreted the mean-field transition 
temperature as the phase transition temperature to a 
superconducting state, as valid in BCS theory. How- 
ever, with an underlying real-space short-range pairing 
and, in many cases, small Fermi surfaces in graphitic 
systems, it is reasonable to question if mean-field the- 
ory is able to estimate the superconducting transition 
temperature and not just to give a temperature scale for 
onset of incoherent mean-field electron pairing. The thus 
formed bound-electron pairs might eventually condense 
but then through a Bose-Einstein condensation (BEC) 
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and then at an often significantly lower temperature. 
The so-called BCS-BEC crossover problem (for reviews 
see e.g. lol"l | 52l | 53l) describes the problem of discerning the 
appropriate theory for a specific system with effective 
fcrmion attraction. Pistolesi et a/. 54 ' 55 showed that the 
parameter kp^, which measures the coherence length £ 
relative to the interparticle distance, is a natural param- 
eter to follow the evolution from BCS to BEC theory. 
They also showed, by using BCS theory, that it can be 
approximated using the relationship T c /Tp rs 0.4/(fcp£), 
where Tp is the Fermi temperature. This relationship 
will work for Dirac fermions as well if identifying kp by 
Ep = hvfkp, where Ep is the Fermi energy and v/ is the 
Fermi-Dirac velocity. Without going beyond mean-field 
theory we therefore estimate = OA/ (T c /Tf) ~ 300 
for the C2S system. Following Ref. [5J, that places it in 
an Uemura plot 55 significantly closer to the classical BCS 
conductors than the high-T c exotic superconductors. For 
the GICs, the steep increase in T c with only moderate 
increase in doping level cause them to fall within criti- 
cal distance of the crossover to BEC theory at kp£ ~ 2ir. 
This estimate strongly questions our mean-field T c as the 
true superconducting temperature for the GICs as well 
as other potential graphitic systems with a high level of 
doping and thus a high predicted T c , In addition, for d- 
wave symmetry the effective electron pair size is always 
finite and it has been shown that in this case the super- 
conducting bosonic regime is in fact never reached i 53 i 57 
In combination the <i-wave symmetry and the closeness 
to a BEC-like theory might strongly influence the value 
of T c for highly doped graphite bands, possibly to the 
extent as to completely suppress superconductivity. In- 
triguingly, however, it should still be possible to observe 
consequences of the effective pairing attraction even in 
these materials. For the GICs, however, this is under 



the assumption that the large metallic s-band does not 
suppress correlations all-together as argued above. 

It might also be worth commenting on the two- 
dimensionality of our model. In graphite there exists a 
non-zero interplanar coupling which should make a mean- 
field transition approximately valid. In the case of inter- 
face superconductivity, such as suggested above for the 
graphite-sulfur composites, the transition should instead 
be of the Kosterlitz-Thouless (KT) type found in thin 
films. Here the mean-field transition temperature will 
give the temperature scale at which the amplitude of the 
order parameter becomes non-zero whereas phase coher- 
ence occurs only at the lower KT transition temperature. 

In summary we have demonstrated the possibility to 
achieve d-wave superconductivity in graphite layers at 
low but non-zero doping levels through an effective favor- 
ing of spin-singlet nearest neighbor bonds. The d-wave 
state is TRS breaking, opening the possibility of coexist- 
ing ferromagnetic and superconducting correlations. We 
have also shown by DFT calculations that sulfur can in- 
duce high enough hole doping levels in graphite to render 
the interface in graphite-sulfur composites superconduct- 
ing through our proposed mechanism. Finally, we have 
noted the possible failure of mean-field theory to esti- 
mate T c for the GICs and other highly doped graphitic 
systems. 
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